Surface states and local spin susceptibility in doped three-dimensional topological 
insulators with odd-parity superconducting pairing symmetry 
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We investigate characteristic features in the spin response of doped three-dimensional topological 
insulators with odd-parity unequal-spin superconducting pairing, which are predicted to have gapless 
Majorana surface modes. These Majorana modes contribute to the local spin susceptibility, giving 
rise to a characteristic temperature behavior of the Knight shift and the spin-lattice relaxation time 
in magnetic resonance experiments. Because of their different decay lengths, the Majorana modes 
can be observed and clearly distinguished from the Dirac modes of the topological insulator by local 
probes which allow for a depth-controlled study of the electron spins on the nanometer length scale. 
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Topological insulators (TIs) are time-reversal invariant 
systems with gapped bulk and protected massless Dirac 
modes at the surface (TJ [2] . Very soon after the discov- 
ery of TIs, theorists generalized this concept to topolog- 
ical superconductors (TSCs) which are characterized by 
a superconducting (SC) gap in the bulk and protected 
gapless Majorana surface modes [TH5]. Semiconductors 
like the bismuth chalcogenides with strong spin-orbit cou- 
pling and a Fermi surface centered at the time-reversal- 
invariant momentum, are of particular interest because of 
their single helical Dirac cone at the surface [7]. Copper 
doped Bi2Se3 is an unconventional SC [BUTT] with non- 
trivial surface states and a band structure similar to that 
of Bi2Se3 but with shifted chemical potential, reduced 
Fermi velocity, and enlarged surface dispersion [9]. 

By now, the surface states in Cu x Bi 2 Se3 have been 
probed by photoemission [S] and point contact spec- 
troscopy [1_21 Q2]. Nuclear magnetic resonance (NMR) 
and quadrupole resonance, as well as the electron (ESR) 
and muon spin resonance (/xSR) are another class of pow- 
erful techniques to investigate the electronic properties 
locally. These probes measure the spin susceptibility 
which is directly connected to the spin structure of the 
SC pairing. In this letter, we study characteristic fea- 
tures in the spin response of odd-parity pairing in doped 
TIs and predict clear signatures for the above resonance 
techniques. 

Fu and Berg E] showed that strong spin-orbit coupled 
bands indeed favor an odd-parity interorbital unequal- 
spin pairing |15[ I16j . To gain insight into its topologi- 
cal non-trivial nature, we map this pairing Hamiltonian 
onto the conduction band which yields an effective time- 
reversal invariant p±ip pairing in three dimensions. Be- 
cause of this topology, there is a pair of Majorana zero- 
energy modes (MZM) located at each surface and pro- 
tected by time-reversal symmetry. Additionally there are 
unconventional Surface Andreev bound states (SABS) 
originating from the band inversion as shown in Ref. 
[17j for a linear k ■ p model. In addition to the linear 



momenta, we here consider quadratic momentum terms 
which determine the energy range of coexistence between 
Dirac modes and unconventional SABS, and which may 
give rise to another species of zero-energy SABS. The 
main motivation for introducing the quadratic terms is 
the possibility to study the competition between the dif- 
ferent surface states and the bulk. The coexistence of the 
MZMs and the SABS originating from the band inversion 
gives rise to two characteristic length scales. The Dirac 
modes decay on the nm scale £ whereas the decay length 
£i for the MZMs is hundreds of nm. Hence, the local spin 
susceptibility shows different characteristic behavior in 
the bulk, at the surface, and within £i into the bulk de- 
spite the spinless nature of the MZM. Conventional bulk 
NMR can distinguish between competing pairing sym- 
metries by the characteristic temperature dependence of 
the Knight shift and the spin-lattice relaxation rate. We 
propose that depth controlled probes [TSJ E] allow to 
clearly determine the pairing symmetry and investigate 
the MZMs. Our work is motivated by Cu x Bi2Se3, how- 
ever, our results are more generally relevant for other 
inversion symmetric materials such as the ternary chalco- 
genides [2U] and the PbTe class [22 [52]. Furthermore, our 
findings for doped TIs are complementary to the super- 
fluid 3 He-i? phase 23, 24 where spin relaxation reflects 
the gapless Majorana nature. 

In the following, we consider doped three-dimensional 
TIs described by the low-energy k ■ p Hamiltonian Hti = 
Ek C^K T /(k)C k [3 El, where 

HtiQl) = m (k)a x + v z k z <Ty + v(k x s y -kyS x )a z - (j, (1) 

with m (k) = to + Bik 2 z + B 2 (k 2 x + h%), C k = (c iks ) jlS , 
Bi,v z > 0, and to < 0. Here, the (<7j) denote 
the Pauli matrices for the spin (orbital) degree of free- 
dom, and the operators c;ks annihilate an electron in 
orbital i, with momentum k, and spin s. The uni- 
tary operator exp(i(f>s z /2), where </> is the azimuthal 
angle (k x ,k y ), transforms Hamiltonian Eq. ([!]) onto 
T~LTi{k x ,§,k z ) which reflects the rotational symmetry 
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around the z-axis. The doped charge density determines 
the chemical potential /i, and the Fermi surface is given 
by the equation fi 2 = mj^kp ) + v 2 z k 2 Fz + v 2 (k 2 Fx + k Fy ). 

Possible pairing terms depend on the specific mech- 
anism and the lattice symmetry. For pairing induced 
by on-site density-density interactions, Fu and Berg [TI] 
showed that the spin-orbit coupled bands favor odd- 
parity interorbital unequal-spin pairing, 

H S c = A ^ ( c ikt c 2-ki + cik4.c 2 -ktj + h.c. (2) 

k 

Here, the parity operator V transforms c\\^ s into c 2 -ks 
and vice versa which yields VHsc'P 1-1 = —Hsc- If 
the pairing is induced by long-range interactions, such as 
the electron-phonon interaction, other unequal-spin pair- 
ing channels are also possible [13 [35] . Additionally, we 
consider (i) even-parity, intraorbital Ai(ci^ci^ + c 2 ^c 2 ^), 
(ii) odd-parity, intraorbital A 2 (ci^c\^ — C2fC 2 i), and (hi) 
even-parity, interorbital A 3 (ci^c 2 ^ — Ci^c 2 -f) pairing. 

To study the effect of spin-orbit coupling on the pair- 
ing symmetry, we project Eq. ^ onto the basis («k,r) 
spanned by the conduction band of the bulk TI Eq. ([!]). 
This yields the effective pairing Hamiltonian 

Ev z k z + ivk m °^ 
„ j, v M a k+ a_ k _ + h.c. (3) 

where £$(k) = v l k l + v2 k 2 and fc2 = k l + k l- Tne 
effective pairing Hamiltonian is exact in first order in 
A//x and yields a gapped bulk excitation spectrum with 
quasiparticle gap 2A_E (k)//i. The operators a^ T and 
c*-k,T transform into each other under time-reversal 
while Vol^.tV^ 1 = ra_ ki _ T [25]. The corresponding 
Bogoliubov-de Gennes equations are in the same univer- 
sality class as the ones for two copies of spinlcss SCs with 
opposite chirality, which is known to be a time-reversal 
invariant TSC with MZMs if the chemical potential lies 
within the conduction band [4,26 . From this analogy, we 
expect to find a pair of MZMs for k — and additionally 
a pair of zero-energy SABS whenever mo(kp) = 0. 

The Hamiltonian Eq. ([I]) confined in z-direction has a 
two-dimensional helical massless Dirac cone at the sur- 
face [TJ. In contrast, we also know that the effective pair- 
ing Hamiltonian Eq. ^ yields a pair of helical MZMs at 
the time-reversal invariant momentum k = [4 . Hence, 
we obtain two species of SABSs originating from the 
band inversion and the p ± ip pairing, respectively. We 
find these states by replacing k z — > —id z and solving 
the corresponding Schrddinger equation with boundary 
condition 03^(2 = 0) = %p{z = 0). For the normal con- 
ducting Hamiltonian Eq. (JIJ, we obtain the dispersion 
Ejy T (k) = vk and the surface wave functions 

1> DlT {z,k,<l>) = ^7Jv~ z (e vz ' v %Qi) a ® {l,iTe^) s , (4) 

where v = (l - ,/! + \B x (m + B 2 k 2 ) /v 2 ) /2B 1 . From 
this expression, we find the decay length £0 — 2B\/v z 



for k 2 < (-to - v 2 z /AB 1 )/B 2 and Co = 2B x /v x /(l - 
y/l + AB 1 {m + B 2 k 2 )/v 2 ) for {-m-v 2 z /AB 1 )/B 2 < k 2 < 
—m/B 2 . For k 2 — > —m/B 2 , £ diverges and the surface 
states become bulk states as shown in Fig. [TJa). 

Similarly as for the Dirac modes, we find the k = 
MZMs by solving the Schrddinger equation. Thus, we 
obtain the Majorana dispersion Ej^(k) ~ kvA(m + 
Bxk 2 Fz )jp? near k = 0. The MZMs decay on a char- 
acteristic length £1 as shown in Fig. [ijb) where we plot 
£1 as function of the chemical potential. For fj, > \m\, 
we find £i(aO = ^F,z/A with the Fermi velocity in z- 
direction Vp tZ . Typical values for £1 are by a factor of 
Eq/A ~ fO 2 larger than £ - In contrast to the Dirac 
modes, the k = MZMs exist for all values of the chem- 
ical potential and enter the bulk for fi 00 only. For 
both species of SABSs, the quadratic terms B\ and B 2 
significantly determine the behavior of the decay lengths. 
For B x = B 2 = 0, we obtain £ = v z/\ m \ and £1 = v z /A 
whereas the decay lengths for Bi,B 2 > strongly depend 
on the Fermi velocity and thus on ji. 




FIG. 1. (Color online) Decay lengths of (a) the Dirac mode of 
the TI and (b) the k = MZM of the odd-parity interorbital 
SC with A = 1 meV, m = -0.3 eV, v n = 4 eVA, and B n = 10 
eVA 2 . The lines are defined in panel (a). 

Depending on the doped charge density, the SC state 
could occur with the chemical potential either in the re- 
gion where the Dirac surface states are separated from 
the bulk conduction band [9] or where only the bulk 
states remain. Since the numerator of Eq. ([3| vanishes 
for the Dirac modes with k z — > —iv/v z and fi = vk, 
they are not gapped by Hgc an d yield a ring struc- 
ture of zero-energy SABSs. Higher order terms which 
couple valence and conduction band do not change the 
character of the surface modes qualitatively. The au- 
thors of ref. [17] showed that a branch of SABSs con- 
nects the MZM and the Dirac-like SABS due to the mir- 
ror helicity of the Hamiltonian which here shows up as 
the sign of the mass TO-o(k^). This mass is negative 
near the bottom of the conduction band and changes 
sign for fj 2 = — mv 2 z jB\. For small chemical poten- 
tials, /1 2 < —mv 2 /Bi 1 —mv 2 /B 2 , the mass term is neg- 
ative, and hence, both the MZM and the branch of 
zero-energy SABS originating from the Dirac mode ex- 
ist. On the other hand, for large chemical potentials, 
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—mv1/B\, —mv 2 /B2 < (J. 2 , the mass term mo(kp) is pos- 
itive and only the k = MZM exists. In the regime of in- 
termediate /i, we distinguish two cases. For —mv 2 /Bi < 
/i 2 < —mv 2 /B 2 , there is a momentum hp such that the 
mass term mo(kp) vanishes. This yields another species 
of zero-energy SABS emerging at k = which now car- 
ries the negative velocity from the band inversion as ex- 
pected from Eq. Q and moves towards the Dirac mode 
with increasing fi. Then, both finite momentum SABS 
meet and gap out for /j? — > —mv 2 /B 2 . In contrast, for 
—mv 2 /B 2 < fi 2 < —mv 2 /Bi, the Dirac SABS disap- 
peared in the bulk and is replaced by the unconventional 
SABS discussed by Hsieh and Fu 17 while m + B\k 2 F z is 
negative |25j . Moreover, with increasing chemical poten- 
tial, this finite-fc zero-energy SABS moves towards k = 
and disappears for /i 2 — > —rav 2 jB\. Thus, we conclude 
that the number of species of zero-energy SABS is even 
for m 2 < /i 2 < —mv 2 /Bi and odd for /j, 2 > —mv 2 /Bi. 
Hsieh and Fu [17] do not find the competition between 
these different sectors since the parameters B\ and B 2 in 
their lattice model are small and thus they only consider 
the small chemical potential regime with the k = MZM 
and one species of finite k modes. 

In the following, we investigate signatures of the uncon- 
ventional SABS in the local quasiparticle density of states 
(LDOS) and the local spin response of the TSC con- 
fined in z-direction. For the numerics we discretize the 
z-direction of Hamiltonian Eq. (JlJ with surfaces at z = 
and z = Na by replacing k z ip — > — i(tp n+ i— ?/> n _i)/2a and 
k 2 ^p — > —(tf) n+ i + V'n-i ~ 2ifj n )/a 2 where z — na. Moti- 
vated by the Bi chalcogenides [2] , we use the parameters 
N = 400, a = 6 A, m = -0.3 eV, v = 2v z = 4 eVA, 
B 2 = 2B 1 = 10 eV A 2 , (i = 0.5 eV, and A = 3 meV. 
Here, B 2 is reduced as compared to Bi2Se3 [2] to guaran- 
tee the existence of the Dirac modes and its separation 
from the conduction band as found for Cua;Bi2Se3 [S]. 
From the above analysis we know that the k — MZMs 
and one finite k zero-energy SABS coexist in this regime. 
For abbreviation, we neglect the quadratic momentum 
terms in our analytical results where they mainly renor- 
malize the Fermi velocity. 

In figure [2]ja) we plot the LDOS which shows two 
qualitatively distinct regions. In the bulk we obtain 
N B {e) = 2^E (k F )e/(TTv 2 v z ^/e 2 ~ A 2 E 2 {k F )/y?) with 
quasiparticle gap 2AE (A;_f)//x and sharp coherence peak, 
while there is a finite midgap LDOS at the surface. De- 
pending on whether the Dirac modes already crossed the 
bulk band, we distinguish between fi 2 > v 2 \m\/B 2 where 
only the MZMs appear with 

A/- S (e, z > 0) * e ^"f ' sin 2 (^°) (5) 

and p? < v 2 \m\/B 2 where both MZM and Dirac-like 
SABS exist with 

AAg(£,Z>0)- M 2 + ; 2 2A2 Sln — ( 6 ) 




(o) z/L (d) 



FIG. 2. (Color online) (a) LDOS, (b) Knight shift, and imag- 
inary part of the transverse spin susceptibility as function of 
(c) temperature and (d) excitation energy with T = 0. All 
quantities are normalized to the normal-conducting bulk. 

for e <C A as shown in Fig. [2ja) . The different surface 
states can be clearly distinguished by their decay lengths 
and energy dependencies. Furthermore, the MZMs oscil- 
late with a period v F ,z which is on the nm scale. 

In figure J2jb) we plot the real part of the integrated lo- 
cal spin susceptibility which is proportional to the Knight 
shift K{z) ~ Re J dz'x zz (z,z';q_\\ — 0). Here, the dis- 
tinction between the bulk and the surface is even clearer 
than for the LDOS. Because of the unequal-spin pairing, 
we find a significantly reduced contribution in the bulk 
for T < A. However, K(T — > 0) does not vanish because 
of the strong spin-orbit coupling [27] . In contrast, we 
find a large shift for T A near the surface which is 
even larger than the bulk shift in the normal state be- 
cause of the large midgap LDOS Eq. The tempera- 
ture dependence of the Knight shift is shown in Fig. [3^ a) 
where the light feature is determined by the bulk and 
the weaker lines at larger K by the surface states. The 
Knight shift from the surface states is spread over a wide 
range due to the exponential decrease in the LDOS and 
shows peak-like subfeatures determined by the LDOS os- 
cillations. The lines with largest shift originate from the 
surface where Dirac and MZMs sum up which gives rise 
to a very large LDOS and therefore a large spin response. 

The imaginary part of the transverse spin susceptibility 
Fig. J2jc) is proportional to the NMR rate l/(T x T) ~ 

ImJJqX h( z > z ;q||) [2S]. Similarly to the Knight shift, 

we clearly distinguish between the bulk and the surface. 
The bulk states give rise to an activation law for T <C 
A and the Hebel-Slichter coherence peak for T — > A. 
In contrast, we find a finite T — 0- value for \/{T\T) 
near the surface where the contributions from the Dirac 
and the Majorana-like SABSs can be distinguished by 
their temperature behavior due to the almost constant 
LDOS of the Dirac modes. This temperature dependence 
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is shown in Fig. [3jb) where the surface contribution is 
again spread with peak-like subfeatures. However, the 
rate directly from the surface [z % £o) is much larger 
than the rate from the MZM only (z » ( ) which allows 
to clearly distinguish the MZM from the Dirac modes. 

In figure [2jd) we show the imaginary part of the dy- 
namical transverse spin susceptibility as function of ex- 
citation energy. The spin excitation spectrum shows a 
very different behavior at the surface and in the bulk. 
For us > 2AE (kp)/ n, there is a continuum of bulk 
spin excitations which is sharply bounded from below. 
In contrast, we find edge-edge spin excitations for lj < 
2AEo(kp)/ fj, at the surface and bulk-edge spin excita- 
tions for lj > AEo(kp)/fx. However, the intensity of 
the edge-edge spin excitations is significantly reduced as 
compared to the bulk-edge spin excitations. 
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FIG. 3. (Color online) Distribution of (a) Knight shift and (b) 
spin-lattice relaxation rate as function of temperature. The 
bright features are determined by the bulk, while the broad 
subfeatures are the surface-state response. All quantities are 
normalized to the normal conducting bulk. 

Finally, we discuss experimental measurement schemes 
of the spin susceptibility. In Refs. 2!) -3Tl it has been 
shown that 77 Se and 209 Bi NMR are suitable methods to 
investigate the bulk of Bi2Se3 where the dipole hyperfine 
coupling dominates due to the p-orbital character of the 
Fermi surface. We conclude that for low temperatures, 
NMR can also study the SC phase to determine the pair- 
ing symmetry. In our analysis, we find significant differ- 
ences in the bulk spin response for the competing pairing 
symmetries |25j . The pairing A2 yields a temperature 
independent Knight shift and a NMR rate with power- 
law T 5 , while A 3 shows a broadened Hebel-Slichter peak 
and an enhanced Knight shift because of a strongly mo- 
mentum dependent gap. Hence, bulk NMR can clearly 
distinguish A2 and A3 from A, while A and Ai are quali- 
tatively similar. The response for Ai is isotropic, whereas 
the A-case yields an anisotropic Knight shift [36]. As 
shown in Fig. [3j a more direct way to distinguish A and 
Ai are the SABS for A which can be observed in pow- 
der samples and thin films as additional signals. The 
integrated intensity of the surface signal is reduced by 
a factor of £i/L compared to the bulk signal where L 
is the thickness of the flakes. However, the surface sig- 
nal shows a characteristic temperature behavior with de- 
creasing spin susceptibility for T — > A. 



As shown in Fig. [2jc) , the SABS yield a spatially de- 
pendent relaxation rate up to v z /A ~ 200 nm into the 
sample which is a characteristic signature for the topolog- 
ically non-trivial nature of the SC state [35] • To investi- 
gate this local magnetism, depth resolved techniques such 
as fjSR Q2] and ,3-NMR Pll[3J] are suitable. In particu- 
lar, they could differentiate the regimes with one, two, or 
three zero-energy SABS, which occur depending on the 
relative strength of the linear and quadratic terms in the 
Hamiltonian as discussed above. For Bi 2 Se3, it has been 
shown that the implanted muons most probably stop in 
the van der Waals gap between quintuple layers [33J . De- 
spite our concentration on NMR and /iSR, our results are 
more general and can also be applied to ESR and surface 
sensitive spin- flip Raman scattering [33] [35] . 

Doped TIs usually have SC shielding fractions of 
50% [H] while the rest is normal conducting. Here, bulk 
NMR could detect the mechanisms which determine this 
fraction. Depending on /z in the normal part of the sam- 
ple, there might be additional signals which can be clearly 
distinguished from the SC ones by the very different tem- 
perature behavior. This key information about the dis- 
tribution of dopants is an important step towards the 
understanding of unconventional SC in TIs. In this pa- 
per, we did not take vortices into account which can be 
distinguished from the SABS by their temperature be- 
havior and the different length scales. 

In conclusion, we have proposed a NMR experi- 
ment to determine the pairing symmetry of doped 
three-dimensional TIs and to observe the MZMs in 
the odd-parity interorbital unequal-spin pairing chan- 
nel. Quadratic momentum terms significantly determine 
the character of the different species of unconventional 
SABSs. We have shown that the surface states of these 
systems have two characteristic length scales originating 
from the band inversion and the odd-parity pairing, and 
typically differ by two orders of magnitude. They con- 
tribute to a local spin susceptibility which can be clearly 
distinguished from the bulk by their characteristic tem- 
perature behavior and the large local density of states. 
Furthermore, we emphasize the usefulness of depth con- 
trolled local probes which directly show the existence of 
unconventional SABSs. 

We acknowledge N. Georgieva, J. Haase, and D. Ry- 
bicki for useful discussions, and BMBF for financial sup- 
port. 
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Complementary to the main part of this paper, we 
here investigate the excitation spectra and the spin re- 
sponse of doped TIs with unequal-spin pairing and show 
that they are very different from the odd-parity interor- 
bital case. We consider the (i) even-parity, intraor- 
bital Ai(c^cii + C21-C24.), (ii) odd-parity, intraorbital 
^(ci^cii — C2-1-C24.), and (hi) even-parity, interorbital 
A 3 (ci t c 2 i - C14.C2-1-) pairing symmetries. 

As in the main part, we study the effect of spin-orbit 
coupling on the pairing symmetry by projecting the pair- 
ing terms onto the conduction band of the doped bulk 
TI. The conduction band is described by the operators 

Ok,r = E S aV'r <7 ( k ) c CT,k, s with 



W(k) 



1 NE + vkTe~ ix ' 2 
V V# - vkre^ 2 



-i(0/2+T7r/4) 
i(c/>/2+TTT/4) 



and e* x = (m (k) + iv z k z )/ -^/m^k) + v^k^. The 
operators ak, T satisfy the condition a^,4>+2-K,k z ),r — 
—ce(k,<f>,k z ),T, i- e - the operators are An periodic un- 
der rotation in momentum space. Because of the he- 
lical band structure, the operators ot(k,4>,k z ),T trans- 
form into ct(k,<t>+K-k z ),T under time reversal. Moreover, 
under parity transformation ct(k,4>,k z ).T transforms into 
Ta (k.<p+Tt-k z ).-T- Hence, it is straight forward to verify 



(7) 



G 



the effective pairing 

Hi ~ A 1 2J a(fc,0,fc 2 ) : Ta(fc,0+7r : *;,),r + h.c, (8) 



k,r 



#2 ~ A 2 ^ a(fc,0.fc z ),ra(fc.0 +7 r.fe z ),r + h.c, (9) 



#3 « As 51 



m (k) 



/t 



-'(fe,0,fe z ),ra(fc,0+7r,fc i ),r + ^- c - (10) 



These equal pseudospin pairing terms do not vanish be- 
cause of the 47r periodicity of the operators Q!(fe,^,fe a ) !T - 
Even if these effective pairing symmetries are very sim- 
ilar, we observe differences in the effect on the Dirac 
surface states and characteristic bulk behavior when the 
mass term m (k) changes sign. In figure [4j we plot the 
energy of the lowest quasiparticle level for the SC con- 
fined in z-direction as function of in-plane momentum k 
and chemical potential /i. 
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FIG. 4. (Color online) Lowest quasiparticle energy in units 
of Aj for the doped TI confined in z-direction and with (a) 
Ai, (b) A2, (c) A3, and (d) A pairing as function of in- plane 
momentum k and chemical potential fj,. Note the different 
color scales. 



Even-parity intraorbital pairing Ai 

For Ai, the effective pairing Eq. ^ is intraband with 
the same amplitude for both pseudospins. Hence, we find 
conventional s-wave behavior for a metal with spin-orbit 
coupling. As shown in figure [4^a), the single-particle ex- 
citation spectra in the bulk and at the surface are both 
fully gapped with gap 2Ai. The gap at the surface orig- 
inates from the fact that the helical Dirac modes have 
only contributions from one orbital and the pairing is 
intraorbital. We find that even if the Dirac modes are 
gapped, they do not hybridize with the bulk states be- 
cause there is no particle-hole mixing between the bulk 



and the surface and the gap at the surface arises from 
the pairing between the helical Dirac modes only. 
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FIG. 5. (Color online) Ai-pafring: (a) LDOS, (b) Knight 
shift, and imaginary part of the transverse spin susceptibility 
as function of (c) temperature and (d) excitation energy. All 
plots are normalized to the normal conducting bulk. 

The LDOS is shown in Fig. p^a) with an ordinary s- 
wave A/s(e) = 2/xi?o(fcF)e/(7™ 2 ^z\/e 2 — A~J) in the bulk. 
Depending on whether the Dirac surface states of the TI 
already crossed the bulk states, we distinguish between 
/i 2 > v 2 \m\/B2 with bulk contributions only and /i 2 < 
v 2 \m\/B2 where 



A/"s(e,z>0) 



fie 



-2z/Co 



A 2 



(11) 



In our numerics, we only consider the case /i 2 > v 2 \m\/ B2 
with Dirac modes. 

As shown in Fig. [5^b) , the Knight shift is significantly 
reduced in the bulk for T < Ai due to the SC gap and has 
a finite T = value determined by the strong spin-orbit 
coupling. The temperature dependence of the Knight 
shift distribution is plotted in Fig. [6ja) where the dark 
feature is the signal from the bulk with a characteristic 
decrease for T —> 0. The line with larger shift originates 
from the Dirac surface states and shows qualitatively the 
same temperature dependence as the bulk shift as ex- 
pected from the gapped LDOS. 

Similarly to the Knight shift, we find conventional s- 
wave behavior for 1/(T\T) in the bulk and at the sur- 
face with an activation law for T <C and the Hebel- 
Slichter coherence peak for T — > Ai as shown in Fig.[5jc) 
and Fig. [6jb) . Again, the surface shows qualitatively the 
same behavior as the bulk but with a much larger rate 
due to the large surface LDOS. The transition between 
the characteristic surface and bulk behaviors occurs in a 
depth £o- 

In figure [5jd) we plot the imaginary part of the trans- 
verse spin susceptibility as function of excitation energy. 
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For to > 2Ai, there is a featureless continuum of bulk 
spin excitations which is sharply bounded from below at 
lj = 2Ai because of the quasiparticle gap. Again we find 
a larger spin susceptibility at the surface as compared to 
the bulk as a consequence of the large LDOS from the 
Dirac modes. 
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FIG. 6. (Color online) Ai-pairing: Distribution of the tem- 
perature dependent (a) Knight shift and (b) NMR rate. All 
plots are normalized to the normal conducting bulk. The 
dark features are determined by the bulk, while the subfea- 
tures show the Dirac surface response. 



Odd-parity intraorbital pairing A2 

For A 2 -pairing, the orbitals pair into singlets with a 
relative minus sign. Because of this relative minus sign, 
the effective pairing term Eq. ([9| is linear in k and van- 
ishes for k x = k y = which gives rise to point nodes in 
the bulk single-particle excitation spectrum and a linear 
dispersion for e <C A 2 . However, for the Dirac modes 
with /i = vk, we find as for the even parity case a gap 
2A 2 . 
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FIG. 7. (Color online) A 2 -pairing: (a) LDOS, (b) Knight 
shift, and imaginary part of the transverse spin susceptibility 
as function of (c) temperature and (d) excitation energy. All 
plots are normalized to the normal conducting bulk. 

In figure [TJa) we plot the LDOS for A2-pairing which 
is finite for all e > due to the linear dispersion and 



shows a cusp at e = A 2 . In the bulk, we find for the 
LDOS the conventional result for SCs with point nodes, 



A 2 TTV^V Z 



e + A-2 



A 5 



(12) 



As before, we distinguish between fi 2 > w 2 |m|/i? 2 with 



bulk contributions only and ji 2 < w 2 |m|/i? 2 with 



Af S {e,z>Q) 



fie 



-2*/£o 



A 2 



(13) 



The different energy dependence of the bulk and the sur- 
face dispersion yields strong evidence for different char- 
acteristic temperature behavior in the bulk and at the 
surface. 

As shown in Figs. [7]Jb) and[8ja), the Knight shift is 
essentially independent of temperature even for T < A 2 . 
Here, the weak subfeatures at smaller shift originate from 
the Dirac surface states and show an s-wave behavior 
with reduced Knight shift below T c in contrast to the 
constant bulk value. 

In figure [f^c) we plot the NMR rate which is much 
larger at the surface than in the bulk. The temperature 
dependence of the rate is shown in Fig. [Sj^b) where we 
can clearly distinguish between bulk and surface contri- 
butions. In stark contrast to the constant Knight shift, 
we find for T — > a power law T 5 in the bulk which 
is characteristic for point nodes and an activation law 
exp(— A 2 /T) from the surface contribution. For T —> A 2 
we obtain a Hebel-Slichter peak at the surface whereas 
the bulk yields just a small coherence peak due to the 



broadened LDOS Eq. (12 1 



The imaginary part of the dynamical transverse spin 
susceptibility is shown in Fig. (7^d) as function of the 
excitation energy. For all e there is a continuum of single- 
particle spin excitations with an amplitude that behaves 
as function of energy like J\fg(ui/2, z) in the bulk and like 
(Af B (u}/2, z) +N s (u>/2, z)) 2 at the surface. For u « 2A 2 , 
there is a peak in the spin spectrum resulting from the 
cusp in the LDOS. 
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FIG. 8. (Color online) A2-pairing: Distribution of the tem- 
perature dependent (a) Knight shift and (b) NMR rate. All 
plots are normalized to the normal conducting bulk. The dark 
feature is determined by the bulk, while the subfeatures show 
the Dirac surface response. 



Even-parity interorbital pairing A3 



The effective pairing Hamiltonian for A3 is shown in 



Eq. (fOl with a fully gapped bulk single-particle excita- 



tion spectrum for mo(k^ ) 7^ and a momentum depen- 
dent gap 2A 3 |mo(fc)|//i. Since mo(fc) changes sign when 
k increases because of the band inversion at k = 0, we 
find a gap closing and reopening as function of [L which 
is intricately related to the transition from the TI phase 
into the band insulator phase. Furthermore, A3 does 
not gap the Dirac modes for which m (k, —iv/v z ) = v 
and the contributions from v and — v cancel. How- 
ever, from the quasiparticle energies we see that the bulk 
spectrum shows nodes whenever mo(fc) changes sign as 
shown in Fig. [4jc) . Similarly to the A case, we here 
find three qualitatively distinct sectors with transitions 
at /j, 2 = —mvl/Bi and /j, 2 — -mu 2 /^. The origin of 
the nodes and distinct sectors is again the band inver- 
sion near k = which changes into the trivial band order 
along the gapless line in Fig.Qc). However, in contrast to 
the odd-parity interorbital case, here, the gapless modes 
are bulk modes. In our numerics, we only consider the 
case [i 2 < —mv 2 /Bi < — mv 2 / B 2l where the bulk is fully 
gapped and the Dirac modes are ungapped. 




FIG. 9. (Color online) A 3 -pairing: (a) LDOS, (b) Knight 
shift, and imaginary part of the transverse spin susceptibility 
as function of (c) temperature and (d) excitation energy. All 
plots are normalized to the normal conducting bulk. 

As shown in Fig. Ma), the LDOS in t he bulk A/ S (e) = 
2^,E (kF)e/(7rv 2 v z y / e 2 — A^m^kp) / fJ- 2 ) with a quasi- 
particle gap 2A3|mo(fci?)|//i. From this expression, we 
see that the gap closes when mrj(fc) vanishes and de- 
creases with increasing /i. As before, we distinguish be- 
tween fi 2 > v 2 \m\/B2 where only the bulk contributes 
and /i 2 < v 2 \m\jBi with a gapless surface LDOS 



Af s (e,z>0) 



fie 



ir£, v 2 



(14) 



itatively very different behavior in the bulk and at the 
surface. In the bulk we find a reduced shift for T < A3 
with a finite K(T —> 0) due to the strong spin-orbit cou- 
pling. In contrast, at the surface we find a larger shift 
for T <C A 3 because of the strongly spin-orbit coupled 
massless surface states. The temperature dependence of 



the Knight shift is shown in Fig. 10 a,) where the strong 
feature is determined by the bulk and the second line by 
the Dirac modes. 

In figure [9^c) we find similarly to the Knight shift a 
qualitatively different behavior for the NMR rate in the 
bulk and at the surface. The gapped bulk gives rise to an 
activation law for T< A3 and a Hebel-Slichter coherence 
peak for T — > A3. In contrast, we find a finite T = 0- 
value for 1/{TxT) at the surface which is characteristic 
for metallic states. The temperature dependence is also 
plotted in Fig. [Toj^b) where the rate from the Dirac states 
is much larger than the rate from the bulk because of the 
large gapless LDOS. 

In figure |9jd) we show the imaginary part of the dy- 
namical transverse spin susceptibility as function of ex- 
citation energy. The excitation spectrum shows a very 
different behavior at the surface and in the bulk. For lj > 
2A3|mo(fci?)|//z, there is a continuum of excitations in the 
bulk which is sharply bounded at lo — 2A3|mo(fcF)|//z. 
In addition, we find low-energy spin excitations at the 
surface and bulk-edge excitations for w > A 3 \m (kF)\/ H 
within a length £0 into the sample. 
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In figure [9^b) we plot the Knight shift K(z) with qual- 



FIG. 10. (Color online) A3-pairing: Distribution of the tem- 
perature dependent (a) Knight shift and (b) NMR rate. All 
plots are normalized to the normal conducting bulk. The 
dark features are determined by the bulk, while the subfea- 
tures show the Dirac surface response. 



Odd-parity interorbital pairing A 

In the main part of this paper we showed that for 
fi > \m\ 3> A the odd-parity interorbital pairing Hamil- 
tonian can be mapped onto a time-reversal invariant 
p ± ip pairing. Furthermore, we found that this pair- 
ing symmetry gives rise to a Kramers pair of MZMs 
at the time-reversal invariant momentum fc = and a 
pair of zero-energy SABSs whenever m (k^) = 0. In 
contrast, the bulk single-particle excitation spectrum is 
always fully gapped with weakly momentum dependent 
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gap Aa/1 — mj^kj?)//^ 2 > 0. The aim of this section is 
to get insight into the quasiparticle spectrum for different 
parameter regimes. 

At the surface we find an interesting spectrum as 
shown in Fig. Qd) where we see a pair of k = MZMs for 
all fi and depending on the chemical potential, we find 
three regimes which can be distinguished by the num- 
ber of additional finite-k zero-energy modes. However, 
in contrast to the k — mode, the finite-fc modes are 
not MZMs in the sense that they are equal superposi- 
tions of electron and hole creation operators such that 
7 = 7^. Instead, the Bogoliubov operators for these 
finite-k modes satisfy 7k = j_ k with a finite electron- 
hole imbalance. For /i 2 < —mv 2 /Bi, we find additionally 
to the k — MZM one pair of finite k zero-energy modes 
located at k = fx/v originating from the band inversion 
of the TI. They are protected by the mirror helicity of 
the Hamiltonian and henceforth not gapped as can be 
seen from the effective pairing Eq. (pi). In the follow- 
ing, we consider the two cases —mv^/B\ ^ —mv 2 /B2 
which are shown in Fig. |11| When tuning through the 
point fi 2 = —mvl/Bi < —mv 2 /B 2 , the sign of the mass 
term mo(fc^ z ) changes which gives rise to another finite-k 
zero-energy mode which carries the negative sign origi- 
nating from the band inversion. This additional mode is 
shown as the dark dispersing feature in Fig. |4|d) . The 
momentum of the third mode increases with increasing 
chemical potential like 



B x v 2 - B 2 v 2 z 



(15) 



and for [j 2 = —mv 2 /B 2 the finite-k zero-energy modes 
fuse and gap out due to the vanishing of the band in- 
version. Hence, we find only the k = MZMs for /i 2 > 
—mv 2 /B 2 > —mv 2 /Bi. For —mv 2 /Bi 3> — mv 2 /B 2 , we 
find a different behavior because the Dirac modes en- 
ter the bulk for fi 2 —> —mv 2 /B 2 . Thus, the finite k 
zero mode looses its Dirac character and moves towards 



zero momentum as described by Eq. (15 1. However, this 
mode is still protected by the mirror helicity and can 
only disappear for /i 2 — > —mv 2 jB\. Thus, the phase 
with three zero-energy modes is absent in this case as 
shown in Fig. 11 a) . For —mv 2 /Bi > —mv 2 /B 2 , both 



processes compete. In figure [TTJb) , we show that in this 
case another finite k zero energy mode emerges while the 
Dirac-like zero-energy mode moves towards zero momen- 
tum. When the finite k zero-energy modes cross, they 
gap out. All these findings are in excellent agreement 
with our projection analysis in the main part of this pa- 
per. 




n 0.2 1.3 



0.5 




(a) 



0.25 0.5 0.75 1.0 1.25 
k/(-m/B,) 1 ' 2 



0.5 



(b) 



0.25 0.5 0.75 1.0 1.25 
k/(-m/BJ 1 ' 2 



FIG. 11. (Color online) Lowest quasiparticle energy in units 
of A with v = 4 eV A and m — —0.3 eV. (a) Regime with 
v 2 /B 2 < vllBx and v z = 2 eV A, B a = 5 eVA 2 , B 2 = 40 
eVA 2 ; (b) Regime with v 2 / B 2 = v\jB x and v z = 4 eV A, 
Bi = B 2 = 10 eVA 2 . 



Knight shift and NMR rate in the bulk 



In figure |l2|(a) and (b), we summarize and compare 
the results for the spin response of the bulk for the above 
introduced pairing symmetries. In this section we assume 
that the order parameter A, is temperature dependent 
with Ai(T) = Aii/l - T/1.76A 4 . 

For the Knight shift Fig. |l~2|a) , we find three different 
behaviors below T c . All pairing terms yield a finite value 
for the zero-temperature spin susceptibility because of 
the strong spin-orbit coupling [57j. However, for A, Ai, 
and A3 we find a decrease of the Knight shift below T c 
with a similar functional behavior for Ai and A which 
differ by K^/i for T = 0. For A3 we observe a qual- 
itatively similar temperature dependence but the func- 
tional form is very different as for A and Ai because of 
the strongly momentum dependent quasiparticle gap. In 
stark contrast to these reduced Knight shifts, we find no 
change of the Knight shift for A2 below T c . This also 
shows that a constant Knight shift for T < T c is not a 
clear signature for triplet pairing symmetry and can also 
appear for singlet pairing symmetries. 
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FIG. 12. (Color online) Temperature dependence of (a) the 
Knight shift and (b) the spin-lattice relaxation rate in the 
bulk for the different pairing symmetries with A; = 1 meV. 
The lines are defined in panel (a). 

In the NMR rate Fig. [l"2{h), the differences between 
the pairing symmetries are not as strong as for the Knight 
shift. For A and Ai, we find conventional s-wave type 
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behavior with an activation law for T — ¥ and a sharp 
Hebel-Slichter coherence peak for T — > T c . However, we 
can again clearly distinguish A 2 and A 3 from A and Ai 
by their functional dependence. As shown above, we find 
a power-law T 5 behavior for A 2 due to the point nodes 
and for A3, the curve is compressed towards T = due 
to the momentum-dependent gap and shows no sharp 
coherence peak. 

Hence, we conclude that it is clearly possible in an 



experiment to distinguish A and Ai from other unequal- 
spin pairing symmetries by conventional bulk NMR ex- 
periments while for the distinction between A and Ai 
quantitative comparisons are necessary. However, an 
even clearer distinction between Ai and A can be ob- 
tained by probing the existence of SABSs in powder sam- 
ples or thin films, or even better depth controlled local 
probes such as /xSR and /3-NMR. 



